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Abstract 

Collective X-ray Thomson scattering has become a versatile tool for the diagnostics of dense 
plasmas. Assuming homogeneous density and temperature throughout the target sample, these 
parameters can be determined directly from the plasmon dispersion and the ratio of plasmon 
amplitudes via detailed balance. In inhomogeneous media, the scattering signal is an average of 
the density and temperature dependent scattering cross-section weighted with the density and 
temperature profiles. We analyse Thomson scattering spectra in the XUV range from near solid 
density hydrogen targets generated by free electron laser radiation. The influence of plasma 
inhomogeneities on the scattering spectrum is investigated by comparing density and temperature 
averaged scattering signals to calculations assuming homogeneous targets. We found discrepancies 
larger than 10% between the mean electron density and the effective density as well as between 
the mean temperature and the effective temperature. 
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1. Introduction 

Recently, X-ray Thomson scattering has 
been demonstrated to be a reliable and versa- 
tile tool for the diagnostics of dense plasmas [l[ . 
Examples for applications cover a large range 
from laser produced plasmas at temperatures of 
several keV and densities between solid density 
and a few percent of solid density (hot dense 
plasmas) 0, H via so-called warm dense mat- 
ter (WDM) US with temperatures of several 
cV and densities close to solid density @, 0] , up 
to compressed matter at temperatures between 
0.1 eV and several 10 eV Hi]. 

The determination of the plasma properties, 
such as the plasma density, the plasma tem- 
perature, and the plasma composition, requires 
the knowledge of the density and temperature 
dependent scattering cross-section. As will be 
outlined in Sec. [H the cross-section for Thom- 
son scattering is usually expressed via the dy- 
namical structure factor (DSF) 5 ee (k, ui) [1 01 ] . 



which is a function of the scattering wavevec- 
tor k and the frequency shift uj, convoluted 
with the density and temperature distribution 
of the target. In this way, plasma inhomo- 
geneities can be taken into account in the sim- 
ulation of the scattering spectrum, which is the 
main purpose of this paper. Such analysis has 
been performed earlier by Baldis et al. ll| for 



the case of VUV photons interacting with solid 
aluminum targets. Here, the scenario of XUV 
photons delivered by the free electron laser fa- 
cility at DESY-Hamburg (FLASH) [HI scat- 
tering on liquid hydrogen at near solid density 
(mass density p = 0.088 g/cm 3 ) as proposed by 
H611 et al. [13], will be considered. 



Theory of Thomson Scattering and 
the dynamic structure factor 
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The scattered power per solid angle dil = 
sin $ di} dip and unit frequency interval du is 
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given by the expression 11 1 
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(Tt = 6.65 x 10 _24 cm 2 is the Thomson cross- 
section, ki and kf are the initial and the fi- 
nal photon wavevector, respectively. By k = 
kj — kf, we denote the scattering wavevector. 
From the conservation of momentum, we ob- 
tain k — 47r sin(i9/2)/Ao, with Ao being the 
probe wavelength. I(r) is the r dependent 
power density of incoming photons, while A rac j 
is the radiated surface of the target. Via the 
convolution of the DSF with an appropriately 
chosen instrumental function G^ w (w), the fi- 
nite resolution of the spectrometer as well as 
the probe's bandwidth is taken into account. 
Here, we use the normalized Gaussian distri- 
bution 



1 



2ttct 



exp 



Aw = 2y2 In 2 o is the full width at half maxi- 
mum (FWHM) of the instrumental function. 

Finally, the plasma inhomogeneities are 
taken into account by averaging the DSF with 
the electron density profile n e (r) and the tem- 
perature profile T(r). 

It is convenient to separate the DSF into 
contributions from free-free, free-bound and 
bound-bound correlations as proposed by Chi- 
hara [10 ] . 

S ee (k,u;) = |/i(fc) + g(fc)| 2 S u (k,w) 

/OO 
dtJ 5 c (k, a;)5' s (k, to — uj 1 ) 
-oo 

+ Z f S e ° e (k,u,) . (3) 

The first part gives the correlation of electrons 
that are weakly and tightly bound to the ions 
and follow the ion's movement adiabatically. 
The amplitude is determined by the atomic 
form factor /;(£;), i.e. the charge distribution of 
the electrons in the valence shell orbitals and 
the screening cloud q{k) which gives the dis- 
tribution of quasi-free electrons screening the 
ion's charge [141 ]. 



The second term contains the contribution of 
core electrons, S c (k,uj) and describes Raman 
type transitions of inner shell electrons to the 
continuum, modulated by the ion's movement 
which is contained in S 8 (k, ui) (H)]. 

Finally, Sg e (k, oS) is the free electron contri- 
bution. It determines the behaviour of the to- 
tal electron structure factor at frequencies close 
to the plasma frequency. Since we focus on 
the application of Thomson scattering to di- 
agnostics of dense plasmas via analysis of the 
plasmon feature, the third term in Eq. (|3|) will 
be discussed in more detail in the following. 
This contribution is calculated in the Born- 
Mermin approximation (l6| . thereby including 
collisions among electrons and ions in second 
order Born approximation. Higher order terms 
can be treated via the t-matrix and inclusion 
of dynamical screening, see Ref. lTf. The in- 



fluence of collisions on the collective scatter- 
ing spectrum has been thoroughly discussed by 
some of the authors in Ref. [18]. In this paper, 
we will concentrate on the impact of the den- 
sity and temperature profiles on the scattering 
signal. 

3. Temperature and density profiles 

We consider the interaction of cryogenic 
(T ~ 20 K) liquid hydrogen droplets of about 
30 fxm diameter with XUV free electron laser 
radiation delivered by FLASH at 13.5 nm wave- 
length with a pulse length of 30 fs and 50 /iJ 
pulse energy. In order to infer the absorp- 
tion characteristics of the target, we calcu- 
late the inverse absorption length a(ui) of cryo- 
genic hydrogen (T = 20 K) via ab-initio quan- 
tum molecular dynamic simulations using the 
Kubo-Greenwood formula [n| • The inverse ab- 
sorption length as a function of the photon en- 
ergy is shown in Fig. [T] At the FLASH en- 
ergy (92 eV, indicated by the red arrow) the 
absorption length is 15 /an, in excellent agree- 
ment with tabulated x-ray absorption data [20( , 
which is in the same order of magnitude as 
the diameter of the hydrogen droplet. Thus, 
FLASH photons can penetrate deeply into the 
target. 

The interaction of FLASH photons with 
the hydrogen target was simulated using the 
radiation-hydrodynamic simulation code HE- 
LIOS [2l[. HELIOS features a Lagrangian ref- 
erence frame, separate ion and electron temper- 
atures, and flux-limited Spitzer thermal con- 
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Figure 1: (Color online) Ab-initio simulation of the in- 
verse absorption length for cryogenic hydrogen at solid 
density as a function of the photon energy. The red 
arrow indicates the FLASH energy (92 eV). 

ductivity. It allows for deposition of laser 
energy via inverse bremsstrahlung as well as 
bound-bound and bound-free transitions, using 
a SESAME-like equation of state. 

Figure [2] shows the free electron density pro- 
file (a) and temperature profile (b) (red solid 
curves). The density profile decreases exponen- 
tially from n c = 2.2 x 10 22 cm~ 3 at the irradi- 
ated surface (x = 0) to n e = 1.0 x 10 20 cm" 3 at 
the rear surface (x — 3.0 /im). The 1/c-dccay 
length of the exponential profile is roughly 
11 /im, which is consistent with the absorption 
length obtained from the QMD simulation, c.f. 
Fig. [TJ i.e. the absorption is mainly due to 
bound-free transitions. 

The plasma temperature decreases from T = 
2eV to T = 1.1 eV. For comparison, the 
dashed lines indicate the mean electron density 
n c = 7.0 x 10 21 cm~ 3 and mean temperature 
T = 1.5 eV, respectively. 

4. Results 

In Fig. [31 we show calculations of the scat- 
tered intensity as a function of the photon en- 
ergy shift with respect to the FLASH photon 
energy, assuming FLASH photons scattering 
from partially ionized hydrogen under an angle 
of d = 90°. The thin curves show calculations 
of the scattering spectrum using a single den- 
sity and temperature combination correspond- 
ing to a particular position x in the HELIOS 
profiles, see Fig. [2] Scattering from plasma re- 
gions close to the front surface (high density, 
high temperature) results in a comparatively 
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Figure 2: (Color online) Density (a) and temperature 
(b) profile from HELIOS simulation of liquid hydrogen 
droplet irradiated by XUV-FEL pulses, Ao = 13.5 nm, 
pulse length 30 fs, pulse energy 50 fiJ (bold red curves). 
The dashed curves indicate the corresponding mean 
density and temperature. 



large plasmon resonance energy, while scatter- 
ing from the rear surface (low density, low tem- 
perature) gives rise to small plasmon energies, 
as indicated by the arrows and the correspond- 
ing electron densities. Finally, the bold red 
curve gives the density and temperature aver- 
aged scattered intensity, c.f. Eq. |T]). Note, 
that since the DSF is weighted with the local 
density n e (r), i.e. with the number of scatter- 
ing electrons in the unit volume d 3 r, as well 
as with the number of photons in that volume, 
regions close to the irradiated surface with in- 
creased density contribute more strongly to the 
spectrum than the dilute parts located towards 
the rear surface. This results in a scattering 
spectrum that is significantly more broadened 
than the original spectra taken at constant den- 
sity and temperature. 
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Figure 3: (Color online) Thin curves: Scattered inten- 
sity spectra at constant density and temperature corre- 
sponding to the HELIOS profiles (Fig. [2}. Bold red 
curve: Density and temperature averaged scattering 
spectrum. Scattering angle •& = 90° . 



We compare the averaged spectrum to a cal- 
culation for the DSF using homogeneous den- 
sity and temperature distributions. The re- 
sult is shown in Fig. |U The bold red curve 
is again the DSF averaged with the HELIOS 
profiles as in Fig. G2 the black dashed curve 
gives the DSF using the effective electron den- 
sity n c ff = 8.3 x 10 20 cm~ 3 and effective tem- 
perature T e ff = 1.7 eV, which are inferred via 
the peak position of the plasmon resonances 
and the detailed balance analysis of the profile- 
averaged spectra (red curve). Both, effective 
density and effective temperature deviate sig- 
nificantly (i.e. more than 10%) from the corre- 
sponding mean values n = 7.0 x 10 20 cm -3 and 
T = 1.5 eV. Also note that in different setups, 
e.g. excitation of the plasma via intense opti- 
cal laser radiation, stronger density and tem- 
perature gradients are observed [22|, and even 
larger deviations between the effective values 
and the mean values for n and T may occur. A 
further interesting observation is that the peaks 
of the averaged spectrum (Fig. |4j are not sym- 
metrically located with respect to the incident 
photon energy as would be the case for a homo- 
geneous temperature and density distribution 
in the target. 



5. Conclusions 

We simulated Thomson scattering spectra 
including density and temperature gradients 
in the target. By comparing these spectra 
to calculations that assume homogeneous den- 



Figure 4: (Color) Density and temperature averaged 
scattering spectrum (bold red curve) compared to a 
calculation using the effective density and temperature 
(dashed curve) obtained from the plasmon peak posi- 
tions and detailed balance analysis from the profile av- 
eraged spectrum. Scattering angle 1? = 90° . 



sity and temperature distributions, we found 
significant differences between these methods. 
In particular, the free electron density that is 
inferred from the simulated spectrum via the 
plasmon dispersion, as well as the effective tem- 
perature, inferred via detailed balance analysis 
of the simulated spectrum may not be inter- 
preted as the mean values of these observables. 
In the scenario studied here (FLASH photons 
at 13.5 nm wavelenght scattering on near solid 
density hydrogen heated by the same photons) , 
both n e g and T e g differ by more than 10% from 
the corresponding mean value. 

Besides being important for the determina- 
tion of density and temperature as demon- 
strated in this paper, deconvolution of the scat- 
tering spectra and the profiles is also crucial 
if one is interested in the broadening of the 
plasmon satellites due to collisions and Landau 
damping 



181 ] . Collisional broadening of plas- 



mons gives valuable information about the elec- 
trical and thermal conductivities in the plasma 
and is thus an important feature of the scat- 
tering spectrum. However, neglecting the ho- 
mogeneous broadening due to the density and 
temperature gradients may lead to significant 
underestimation of the conductivity. Thus, 
gaining information about the density and tem- 
perature profile of the target, e.g. via radiation 
hydrodynamic simulations, is indispensable for 
the application of Thomson scattering spectra 
to plasma diagnostics. 
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